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Abstract 

Very Long Baseline Interferometry (VLBI) Intensives are one-hour single baseline sessions to provide 
Universal Time (UT1) in near real-time up to a delay of three days if a site is not e-transferring the 
observational data. Due to the importance of UT1 estimates for the prediction of Earth orientation 
parameters, as well as any kind of navigation on Earth or in space, there is not only the need to improve 
the timeliness of the results but also their accuracy. We identify the asymmetry of the tropospheric 
delays as the major error source, and we provide two strategies to improve the results, in particular of 
those Intensives which include the station Tsukuba in Japan with its large tropospheric variation. We 
find an improvement when (1) using ray-traced delays from a numerical weather model, and (2) when 
estimating tropospheric gradients within the analysis of Intensive sessions. The improvement is shown 
in terms of reduction of mis of length-of-day estimates w.r.t. those derived from Global Positioning 
System observations. 


1. Introduction 

Very Long Baseline Interferometry (VLBI) is the primary technique for the determination of 
the Earth Orientation Parameters (EOP), in particular for nutation and Universal Time (UT1). 
Typically, 24-hour VLBI sessions with five to eight participating stations are observed on about 
three days per week to provide the EOP with a latency of about two weeks and a UT1 accuracy 
of approximately 6 to 7 /us. Additionally, so-called one-hour Intensive sessions are observed every 
day to provide UT1 with latencies between 3 minutes and 3 days. These Intensive sessions usually 
only contain single baseline observations and yield UT1 estimates with an accuracy of about 15 /is 
when compared to the IERS 05 C04 series (Bizouard and Gambis, 2009, [1]) from the International 
Earth Rotation and Reference Systems Service. 

So-called INTI sessions are observed from Monday to Friday on the baseline Wettzell (Ger- 
many) - Kokee Park (Hawaii, U.S.A.) at 18:30 UT, INT2 sessions from Saturday to Sunday on 
the baseline Wettzell - Tsukuba (Japan) at 7:30 UT, and additionally INT3 sessions on Monday 
at 7:00 UT between Wettzell, Tsukuba, and Ny-Alesund (Spitsbergen, Norway). UT1 estimates 
in near-real time are crucial for the prediction of EOP, and EOP in real-time are needed for any 
kind of navigation and positioning on Earth and in space. In particular, the prediction of Global 
Navigation Satellite Systems (GNSS) orbits benefit from better EOP prediction. Luzum and Noth- 
nagel (2010, [8]) have shown that there is an improvement of 50% in rapid combination and 20% 
in short term prediction if UT1 estimates are available in near real-time, which is the case for 
INT3 sessions and in the foreseeable future also for INTI and INT2 sessions. In this paper, we do 
not assess questions of timeliness, but we discuss the accuracy of UT1 estimates from Intensive 
sessions. 
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2. VLBI Analysis 

We analyzed all Intensive sessions of the International VLBI Service for Geodesy and As- 
trometry (IVS; Schliiter and Behrend, 2007, [12]) from 2007.0 to 2010.0 with the Vienna VLBI 
Software (VieVS, Bohm et al., 2010a, [3]). Station and source coordinates were fixed to VTRF2008 
(Bockmann et al., 2010, [2]) and ICRF2 (Fey et al., 2009, [6]), respectively. We fixed nutation 
offsets and polar motion to the values provided in the IERS 05 C04 series and modeled the high- 
frequency variations in polar motion as recommended by the IERS Conventions 2003 (McCarthy 
and Petit, 2004, [9]). Five parameters were estimated for each single baseline Intensive session: 
offset and rate between the clocks, one zenith wet delay offset at each station, and one UT1 offset 
parameter. In the standard solution, no gradients were applied (neither a priori nor as estimated 
parameters), nor did we use a cutoff elevation angle or downweighting of low observations. Figure 
1 shows UT1 estimates for the INTI and INT2/3 sessions. (We treat the INT2 and INT3 sessions 
as one common series INT2/3). 


Figure 1. UT1 estimates and formal errors from IVS Intensive sessions from 2007.0 to 2010.0 w.r.t. IERS 
05 C04. 

We analyzed 675 INTI sessions and 393 INT2/3 sessions from January 2007 until December 
2009, from which we had to remove 17 and 9 sessions as outliers, respectively (with estimates larger 
than ± 50 /x s or uncertainties larger than 100 /x s). The median biases are 2.4 and 3.7 fxs, and the 
rms values of the UT1 estimates are 15.0 and 14.7 /xs w.r.t. the IERS 05 C04 series. On the other 
hand, the formal uncertainties are 10.9 and 7.3 /xs, which suggests that modeling deficiencies exist 
and need to be investigated. 

A discussion of the influence of nutation and polar motion errors on UT1 estimates from 
Intensive sessions was provided by Nothnagel and Schnell (2008, [11]) who report maximum values 
of 30 /xs in UT1 per mas in nutation or polar motion. However, when using sophisticated forecasts 
of daily polar motion (e.g., from the International GNSS Service IGS) and nutation offsets (e.g., 
from the IERS), this is a minor error source. Presently more important is the high-frequency model 
for polar motion, which has a significant impact on UT1 estimates from Intensive sessions. Figure 2 
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shows the differences in UT1 when using the empirical model as provided by Englich et al. (2008, 
[5]) w.r.t. to applying the model recommended by the IERS 2003 Conventions (McCarthy and 
Petit, 2004, [9]). In recent years, general awareness has grown that there is a need to replace the 
IERS 2003 model (‘Eanes model’) by an improved high-frequency model which is either determined 
empirically from observations, from improved ocean models, or from a combination of both. With 
such a new model, the a priori high-frequency polar motion will contribute less as an error source 
to the determination of UT1 from Intensive sessions. 



Figure 2. UT1 differences in /is when using the empirical model by Englich et al. (2003) for polar motion 
instead of the IERS 2003 ‘Eanes’ model. 

Unlike the error sources described above, asymmetries in the tropospheric delays cannot be 
accounted for easily. Bohm et al. (2010b, [4]) show that the influence of those asymmetries of the 
tropospheric delays on UT1 estimates from Intensive sessions is about ±10 fis but can be as large 
as 50 /j,s during extreme weather conditions. As a rule of thumb, an unmodeled east gradient of 
1 mm (= sum of total east gradients at both stations) causes UT1 to change by about 15 /rs. 
Considering that gradients can reach values of up to 3 mm (corresponding to about 0.3 m delay 
at 5 degrees elevation) this is a very large error source, which needs to be considered. Bohm et al. 
(2010b, [4]) investigated the application of external information about tropospheric delays for the 
analysis of INT2 sessions, and they found a slight improvement when applying ray-traced delays 
at Tsukuba (Hobiger et al., 2008, [7]). Since the IERS 05 C04 series is not good for comparison 
as they already contain information from a standard analysis of Intensive sessions (Bizouard and 
Gambis, 2009, [1]), Bohm et al. (2010b, [4]) compared length-of-day (LOD) between two Intensive 
sessions on consecutive days to LOD derived from Global Positioning System (GPS) observations, 
and the rms of LOD decreased from 23.7 to 22.8 jus for a sample of 70 LOD values when using 
ray-traced delays. 
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3. Estimation of Gradients 

Nilsson et al. (2010, [10]) investigated single baseline observations over time spans of two hours 
that were extracted from CONT08, which was a special VLBI campaign over 15 days in August 
2008 to demonstrate the best accuracy that can presently be achieved with VLBI. This procedure 
allowed Nilsson et al. (2010, [10]) not only to have optimum reference values for UT1 from the 
full CONT08 solution but also to have good a priori values for zenith wet delays or gradients. 
However, one of their findings was that there is an improvement in the estimation of UT1 when 
gradients are estimated from those single baseline two-hour sessions, in particular when Tsukuba 
is one of the observing stations. 

Although the number of observations is small (below 30), we followed the findings by Nilsson et 
al. (2010, [10]) and tested the estimation of gradients in ‘real’ IVS Intensive sessions. In a first run, 
we estimated north and east gradient offsets which are constrained to zero by ±1 mm. In a second 
test, we only estimated east gradient offsets (constrained to zero by ±1 nun), because the east 
component is most important for UT1. The constraints are necessary to avoid singularity of the 
normal equation matrix. As already mentioned above, a quality assessment of the UT1 estimates 
cannot be done by comparison to the IERS 05 C04 values because the latter already contain UT1 
information from a standard solution of IVS Intensive sessions without the estimation of gradients. 
Thus, we again compare LOD values from Intensive sessions on consecutive days to LOD values 
from the GPS solution as provided by the Center for Orbit Determination (CODE) in Bern. Table 
1 shows that we find a considerable improvement when estimating gradients in the INT2 sessions, 
and even more improvement when only estimating east gradients. The reduction in mis of LOD 
values is from 24.1 to 21.6 n s when estimating east gradients at Tsukuba and Wettzell in the INT2 
sessions. On the other hand, there is hardly any improvement when estimating gradients in INTI 
sessions. This is due to the fact that the largest gradients are at Tsukuba, which also confirms the 
findings by Nilsson et al. (2010, [10]). 


Table 1. Rms of UT1 and LOD w.r.t. IERS 05 C04 and GPS, respectively. Considerable improvement is 
found for the rms of LOD w.r.t. LOD from GPS when only east gradients are estimated. 



rms w.r.t. UT1 from IERS 05 C04 

rms w.r.t. LOD from GPS 


INTI 

INT2/3 

INTI 

INT2/3 

samples 

657 

384 

258 

224 

no gradients 

15.0 //s 

14.7 /us 

27.6 fis 

24.i n s 

north and east gradients 

16.1 /us 

15.6 /us 

27.9 ns 

22.3 /us 

east gradients 

16.2 /us 

15.9 /us 

27.0 n s 

21.6 /us 


4. Summary and Outlook 

We discussed possible error sources for the estimation of Universal Time (UT1) from Inten- 
sive sessions, and we identified the asymmetry of tropospheric delays as one — if not the most 
important — of these error sources. Although the number of observations is typically small (be- 
low 30) we found an improvement when estimating (constrained east) gradients in INT2 sessions 
including Tsukuba. In future investigations we will determine the improvement when estimating 


254 


IVS 2010 General Meeting Proceedings 




Johannes Bohm et al.: UT1 from VLBI Intensives 


east gradients only at one station (e.g., Tsukuba). Complementarily, we recommend to continue 
the research in using ray-traced delays for the analysis of IVS Intensive sessions. In future accuracy 
assessments we plan to use different UT1 reference series, e.g., a Kalman Filter series which does 
not include information from IVS Intensive sessions (Ray, 2010, personal communication). 

5. Acknowledgements 

This study made use of VLBI data provided by the International VLBI Service for Geodesy 
and Astrometry (IVS; Schliiter and Behrend, 2007), and the authors would like to thank all IVS 
components. We also acknowledge the Deutsche Forschungsgemeinschaft (DFG) for supporting 
this work (Project SCHUH 1103/3-2). 

References 

[1] Bizouard C, Gambis D (2009) The combined solution C04 for Earth Orientation Parameters consistent 
with International Terrestrial Reference Frame. In Drewes H (ed.) Geodetic Reference Frames, IAG 
Symposium Munich, Germany, 9-14 October 2006, Vol 134, pp 265-270. 

[2] Bockmann S, Artz T, Nothnagel A, (2010) VLBI terrestrial reference frame contributions to ITRF2008, 
J. of Geodesy, 84(3), 201-219, doi:10.1007/s00190-009-0357-7. 

[3] Bohm J, Bohm S, Nilsson T, Pany A, Plank L, Spicakova H, Teke K, Schuh H (2010a) The new Vienna 
VLBI Software. IAG Scientific Assembly Proceedings, Buenos Aires. 

[4] J. Bohm, T. Hobiger, R. Ichikawa, T. Kondo, Y. Koyama, A. Pany, H. Schuh, K. Teke (2010b), 
Asymmetric tropospheric delays from numerical weather models for UT1 determination from VLBI 
Intensive sessions on the baseline Wettzell-Tsukuba, J Geod, 84:319-325, doi:10.1007/s00190-10-0370- 

x. 

[5] Englich S, Heinkelmann R, Schuh H (2008) Re-assessment of Ocean Tidal Terms in High-Frequency 
Earth Rotation Variations Observed by VLBI. The 5th IVS General Meeting Proceedings (ed. by 
Finkelstein A, Behrend D), 2008, pp 314-318. 

[6] Fey A, Gordon D, Jacobs CS (eds.) (2009) The Second Realization of the International Celestial 
Reference Frame by Very Long Baseline Interferometry. Presented on behalf of the IERS / IVS Working 
Group. IERS Technical Note 35, Frankfurt am Main: Verlag des Bundesamts fiir Kartographie und 
Geodasie, 204 p., ISBN 3-89888-918-6 (print version). 

[7] Hobiger T, Ichikawa R, Koyama Y, Kondo T (2008) Fast and accurate ray-tracing algorithms for 
real-time space geodetic applications using numerical weather models. J Geophys Res 113(D20302). 
doi: 10. 1029 /2008 JD010503. 

[8] Luzurn B, Nothnagel A (2010) Improved UT1 predictions through low-latency VLBI observations. J 
Geod. doi:10.1007/s00190-010-0372/8. 

[9] McCarthy D, Petit G (2004) IERS Conventions 2003. IERS Technical Note No. 32, Verlag des Bunde- 
sarntes fiir Kartographie und Geodasie, Frankfurt am Main. 

[10] Nilsson T., Bohm J., Schuh H (2010), Universal Time from VLBI Single-Baseline Observations during 
CONT08, to be submitted to Journal of Geodesy. 

[11] Nothnagel A, Schnell D (2008) The impact of errors in polar motion and nutation on UT1 determina- 
tions from VLBI Intensive observations. J Geod 82(12):863-869. doi:10.1007/s00190-008-0212-2. 

[12] Schliiter W, Behrend D (2007) The International VLBI Service for Geodesy and Astrometry (IVS): 
current capabilities and future prospects. J Geod 81(6-8):379-388. doi:10.1007/s00190-006-0131-z. 


IVS 2010 General Meeting Proceedings 


255 



